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Abstract 

We investigate numerically how accurately one could constrain the strengths of dif- 
ferent short-range contributions to neutrino-less double beta decay in effective field 
theory. Depending on the outcome of near-future experiments yielding information 
on the neutrino masses, the corresponding bounds or estimates can be stronger or 
weaker. A particularly interesting case, resulting in strong bounds, would be a posi- 
tive signal of neutrino-less double beta decay that is consistent with complementary 
information from neutrino oscillation experiments, kinematical determinations of the 
neutrino mass, and measurements of the sum of light neutrino masses from cos- 
mological observations. The keys to more robust bounds are improvements of the 
knowledge of the nuclear physics involved and a better experimental accuracy. 
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1 Introduction 



Within the last decades, it has become more and more apparent that neutrinos can be 
used as one of the major windows to look for physics beyond the Standard Model (SM) of 
elementary particles, which is often dubbed "new physics" . Not only has it been established 
by neutrino oscillation experiments on solar [H [21 El H], atmospheric [5], and artificially 
produced neutrinos [6l [71 Ej that neutrinos are massive, but also the corresponding leptonic 
mixing differs significantly from the quark sector [5]. Apart from their values, one of the 
major questions about neutrino masses is their origin. A particularly attractive possibility 
to explain their smallness is the so-called "seesaw mechanism" [lOl [Til [IS [131 HI]. This 
mechanism, however, intrinsically involves the breaking of lepton number by two units. If 
this breaking was indeed realized in Nature, a process that would have to appear is neutrino- 
less double beta decay (Oi^/3/3), where an even-even nucleus undergoes a transition {A, Z) 
{A, Z + 2) + 2e~, thereby violating lepton number. This process has been intensively 
searched for, but up to now there has not yet been an unambiguous detection [15] (cf., 
however, the claim by part of the Heidelberg- Moscow collaboration in Ref. [T6]). 

Although the "standard" mechanism for Oz//3/9 is the exchange of light Majorana neu- 
trinos [T71 [IHl [19], other mechanisms could very well appear in certain models beyond the 
SM, such as supersymmetry (see, e.g., Refs. [201 [21]) or left-right symmetric models (see, 
e.g., Refs. [221 [23])- To disentangle the various possibilities, it will probably be necessary 
to detect the process for several different isotopes [2^ [25]. 

In this paper, we will investigate point-like contributions of new physics to Oi//3/3, previ- 
ously studied in Refs. [211 [2Z]- Over the last years, there has been considerable improvement 
in the knowledge of the neutrino oscillation parameters [28], and new experiments on neu- 
trino masses and mixings are currently starting or will start in the near future. One of the 
most prominent near future experiments on neutrino-less double beta decay is certainly 
the GERDA experiment [29], which we will take as an example. This experiment will be 
complemented by others aiming to detect the kinematical electron-neutrino effective mass 
(e.g. KATRIN 00]), the cosmological sum S of neutrino masses (e.g. Planck [21]), 
or the leptonic mixing angle ^13 (e.g. Double Chooz [32] )• We therefore feel that it is 
necessary to investigate the impact of these experiments on new physics contributions to 
Oufi(3. Naturally, once this process is observed, there will soon be analyses of the standard 
picture. Later on, however, the question will arise what else we could learn from a positive 
or negative signal of Oi^/3/3. Here, we are going to give an answer to this question in what 
regards the short-range contributions to the process. 

The paper is organized as follows: In Sees. [21 and [31 we will shortly review the nature 
of the short-range contributions to Oiy(3(3, as well as the analysis of a signal under the as- 
sumption of the standard mechanism of light neutrino exchange. Our actual investigations 
will be presented in Sees. [H and [5l where we first define certain benchmark scenarios and 
use them later on to illustrate how strong the bounds on the effective operators of the new 
physics contributions to OufSP could be for certain realistic situations after the future ex- 
periments have yielded results. Depending on if the signal of Ou^f] is positive or negative, 
and also depending on if it is consistent or inconsistent with results from complementary 
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experiments, the bounds can be stronger or weaker. In general, it is hard to derive strict 
statistical bounds, due to the intrinsic difficulties arising from the nuclear physics involved. 
Finally, we present our conclusions in Sec. |6l 

2 Point-like contributions to Oi>'/3/3 

The key point of our analysis is that one can parametrize many different contributions 
of new physics to Oi^/3/3 in a nearly model-independent way by using effective field the- 
ory [331 EH |35]. This is a very powerful approach as long as the detailed spectrum of 
new heavy particles does not play a significant rolell] Depending on if one only modi- 
fies the vertices that couple to the light neutrinos or if the whole decay mechanism is 
altered, one obtains long-range |37] and short-range |26] contributions, respectively. Here, 
we restrict ourselves to the latter ones, since modifications of the vertices only are even 
discussed in some textbooks, see e.g. Ref. [19], and furthermore such modified weak in- 
teraction charged current vertices can be probed much more easily in, e.g., neutron beta 
decay experiments [SH] or in muon conversion experiments 

When looking at the short range contributions, one has to study the following general 
Lagrangian [26] : 

C = ^m;' {eiJJj + 62 J^V^,j + 63 J'V^j + e^J^J^^f + e,,J''Jj^) + h.c. (1) 

Here, Gi;' is the Fermi constant, which is a reflection of the double weak interaction vertex 
involved in the standard light neutrino exchange diagram for Oz//3/3, and is the proton 
mass, which is used to give the Lagrangian the correct mass dimension. In principle, any 
high energy scale could have served this purpose, but a scale of the order of the proton mass 
is already considerably higher than the typical internal momentum in the nucleus, which 
is of the order of 100 MeV [TSj. The strengths of the different operators are parametrized 
by the (generally complex) dimensionless coefficients 6j. The hadronic currents in Eq. ([1]) 
are given by 

JL,R = u{lTl,)d, Jln = ur{lTl,)d, C« = 4b^71(lT75)ci, (2) 
and the leptonic ones are 

3l,r = e{lTl^)e' = 2 e]^(eR,i)^ = (1 T 75) = 2 e^Iil^ien^LT ■ (3) 

For all currents, different chirality structures are permitted. Note that we have also assigned 
the correct chirality to the electron operators in Eq. (|3]), since we will need this for our 
argumentation later on. 

Actually, we could have added some more Lorentz-invariant terms to Eq. ([T]), 

C = ^m;i {e,J^rj,u + e.jr'^j^, + 63 J^^r^J^) , (4) 
"'^See, e.g., Ref. |36| for a case where the detailed heavy spectrum does play a role. 
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where the leptonic tensor current is given by j'^j^ = e| [7^^, 7*^] (1 =F 75) e*^. These terms had 
been included in the Lagrangian in Ref. [26], but were neglected in the final analysis, since 
these authors worked in the s-wave approximation where such contributions vanish. In 
Ref. |27], however, it was pointed out that all operators proportional to e'j^e'^, e| [7^, •ju] e.^, 
and 675! [7^, 7,y] e'^ vanish identically, since the electron fields are Grassmann numbers. 
Therefore, the terms in Eq. (jl]) are not relevant for Oz//3/3. 

The purpose of the effective field theory approach is to provide an easy and unified 
description of new physics. In particular, since there is only a finite set of operators, it 
allows to efficiently compare different underlying models in what regards their interaction 
strengths. Very recent works covering different explicit Oz//3/3 mechanisms and in parti- 
cular their interplay with the standard light neutrino exchange |10l HI] appeared nearly 
simultaneously to this paper. All the mechanisms treated there can be translated into our 
effective language. 

To give one explicit and illustrative example, we take the gluino exchange mechanism 
in _R-parity violating supersymmetry (RPV-SUSY) [101 SH 112] • The goal is to reproduce 
Eq. ([1]) by integrating out all heavy particles, and it is easy to see that the corresponding 
effective Lagrangian is given by 

Cnpv = ^m;i(ef^^)RPv JrJrJr + h.c, (5) 
whose coefficient can be expressed in terms of model parameters as [40] 

1 + (^) , (6) 

where as = gl/ (47r) is the strong fine structure constant, m-^^^^ - are the masses of the up- 
squark, of the down-squark, and of the gluino, respectively, and the parameter A'^^^^ encodes 
the strength of the RPV-SUSY operators. Obviously, if we have any kind of bound on the 
general coefficient ef^'^, this can be translated into a bound on the parameter A'^, and 
hence on the size of the new physics contribution. 

Later on, we will derive bounds on A'^^^^ for all scenarios under consideration, in order to 
illustrate how one can translate our results into bounds on concrete models of new physics. 
In order to arrive at numerical values for A'^^^^^, we will assume that rriuj^ = m^^ = 800 GeV 
and nig = 1000 GeV. These values are good examples for the Minimal Supersymmetric 
Standard Model squark and gluino masses at low energies, when a minimal Supergrav- 
ity scenario is considered to be present at the GUT scale [13]. A generalization of this 
procedure to underlying models different from the RPV-SUSY example is straightforward. 

In Sec. [5|, after having introduced the standard mechanism, we will discuss the actual 
formula for the corresponding decay rate. 



^RRR-. ™« (^111)^ 

= 1 JRPV 



6 G%m^- nig 
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3 The standard analysis of a positive signal 



The most interesting case to study is the one where the GERDA experiment yields a 
positive signal. This corresponds to a measured decay rate Fobs that is different from zero 
and has a standard deviation (XobsS 

Fobs ± (Tabs @ 68 % C.L. 

Assuming that we interpret this measurement as experimental proof of lepton number 
violation, we would, as a next step, try to describe it by assuming the standard light 
neutrino exchange for the mediation of the decay [TS]. In this case, our expectation for the 
decay rate would be 

= Veer, (7) 

where G is a known phase space factor, Ai^'^ is the nominal nuclear matrix element (NME), 
and |mee| is the effective neutrino mass. The latter is a function of the smallest neutrino 
mass m, the mixing angles 612 and ^13, the mass-squared differences Am^i and |Am|^|, as 
well as the two Majorana phases a and (3^ 

To determine the goodness of fit of the measured data with the standard neutrino- 
exchange mechanism, we will construct a x^-function following the procedure outlined in 
Ref. IB], extended by information on m^. In that paper, the authors have used the method 
of the covariance matrix given by 



Sab - ^at-O-Q^exp + 



2 , y^dlam^2 

dxi dxi 



where a,b E {1,2,3}, (7a,exp are the experimental errors of the measured observables, 
Xi G {5^2, -5^3, Am^i, |Am|^|} are the relevant neutrino oscillation parameters, and cTj G 
{^(■^12)5 '^(■^la)' '^(^"^21)5 ^(l^^^iil)} the corresponding standard deviations. Exphc- 
itly, the T^'s are given by 



where 



Ti = ^\mee\, T2 = rua, T3 = S, (9) 



I I I 2 2 I 22 210 I 2 2i/3i 

Irrteel = |miC]^2Ci3 + "^2S]^2Cl3e + m3S;^3e '^1, 

2 222, 222, 22 ^in\ 

= "^lCi2Ci3 + '^2Sl2Ci3 + "^3-5l3' (10) 

E = mi+m2 + m3. 



^Note that, depending on the size of (Jobs, it is possible that this rate is also consistent with a zero decay 
rate at x • cr if Fobs — x ■ cTobs = 0. We will ignore this subtlety here and assume that a positive experimental 
signal means that lepton number violation is experimentally established. This is a reasonable assumption 
in the likely case that we trust the experimentalists conducting the experiment. Furthermore, the analysis 
would become much more subtle when including such pathologic cases, without any significant gain. 

•^Note that \mee\ does not, contrary to some statements in the literature, depend on the Dirac CP- 
phase S, since this phase can be rotated away by redefining the third neutrino mass eigenstate 1/3, see 
Refs. \MM. 
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The above observables are the main sources of information about the neutrino mass that 
we currently have: 

• The effective neutrino mass |mee| is, as just explained, measured in experiments that 
search for Oz//3/3. We will take the GERDA experiment [29] as an example. 

• The effective mass square of the electron neutrino, m|, is measured in kinematical 
experiments on single /9-decay, where KATRIN [30] is currently the most prominent 
upcoming example. 

• The effective sum S of neutrino masses as measured in cosmological observations. 
This can be done through measurements of the cosmic microwave background (CMB) 
radiation, where we take the Planck satellite [31] as example. However, the CMB 
measurements need to be combined with results from other observations such as high 
redshift galaxy surveys, baryon acoustic oscillations, or of Type la supernovae [ITJ HH] 
to yield robust results. 

Using these observables, the full x^-function is given by 

F, a, /3) = min v'^S~^v^ (11) 
Ce[i/VF,v^] 

where v'^ = (Ti - T{'^',T2 - Tf'',T3 - T^^'), and T°''' are the experimental values of the 
measured quantities. In particular, in the case of Ouf^P, we have 

mobs I I /robs 1 / 1 I \ ^obs i i ^ l-'^truel /ir>\ 

Ti = Imeelobs = y ^ |^o^| , ^l,cxp = (T{\mee\) = - | | obs , ^ = |^o^| • (12) 

Here, lA^^^^I is the nominal value of the NME used in the analysis to extract the experimen- 
tal value of the effective mass, whereas lA^trud is its (unknown) true value. Furthermore, 
F parametrizes the uncertainty in the NME, since it forces the x^-minimization to treat all 
values of the NME between \Ai^'^\/\/F and \Ai^'^\y/F on the same footing, corresponding 
to a flat prior. As representative values we will choose F = 1, 2, and 3, where F = 1 
corresponds to a perfect knowledge of the NME (cf. Ref. [S]). Note that, depending on 
which information is available, it might also be required to reduce Eq. fITT]) to incorporate, 
e.g., |mee| and S only, but not m^. From Eq. f lTT]) . it is possible to derive la-ranges for 
the smallest neutrino mass m. 

4 The benchmark scenarios 

The experimental situation in a few years from now is depicted in the left panel of Fig. [H 
which shows the range of the effective mass |mee| in Oz//3/3 as a function of the smallest 
neutrino mass m for both orderings, normal (mi < m2 < m^) and inverted {m^ < mi < 
777,2). The inner (outer) bands show the respective la region for a value of sin^ 613 = (for 
the Double Chooz limit in case of a negative signal, sin^ ^13 < 0.0076 [IH]), while for the 
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'"1.3 [eV] mi_3 [eV] 



Figure 1: The experimental situation in a few years from now and the benchmark scenarios 
(A,B,C) to be considered. 

other oscillation parameters we have used the best-fit values and (symmetrized) standard 
deviations from Ref. [2B]- Obviously, the effect of a tiny but non-zero mixing angle ^13, is 
very small, in agreement with Refs. [501 EI] • fact, for inverted ordering, the outer band 
(lighter yellow) is barely visible. The only effect of a non-zero 6*13 shows up for the region 
of normal hierarchy (normal ordering with very small mi), which cannot be probed in the 
near future. Accordingly, the dependence of the effective mass on the Majorana phase /3 
is very weak [15] . 

We have also depicted the expected sensitivities of future experiments. First of all, the 
GERDA experiment will set different limits on the half-life of Ou^P in three different phases 
I, II, and III, which are given by 3-10^^ years [52], 1.5-10^^ years [52], and 2.8-10^^ years [53] . 
To translate these values into sensitivities on |mee|, one has, in principle, to know the 
correct NME. Since we do not want to enter a discussion about better and worse values of 
this quantity, we have taken, for simplicity, the most optimistic values from Ref. [21] for 
phases I and II. As the reach of phase III is not completely clear yet (and hence, there is, 
to our knowledge, no official number from the GERDA collaboration), we have taken the 
benchmark value |mee| = 0.01 eV, which is in agreement with recent estimates [53] . 

We have also depicted limits from complementary experiments that can yield informa- 
tion on the absolute neutrino mass scale. The only model-independent way to measure 
the neutrino mass is from simple kinematics, which will be performed by the KATRIN 
experiment [30] . Depending on the statistics used, one obtains different predictions for the 
sensitivity of the KATRIN experiment [51], and we have used the slightly more optimistic 
value obtained by using Bayesian statistics, = 0.17 eV (which translates into just the 
same number for the smallest neutrino mass). The Planck experiment is supposed to probe 
the effective sum S of light neutrino masses down to about 0.1 eV [17], which we have taken 
as a benchmark value (and which translates into a limit on the smallest neutrino mass of 
roughly 0.02 eV)0 

^Note that one might have to be careful when taking the cosmological observations at face value [55] . 
since unknown systematic errors might be involved. 
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Based on these numbers, we have chosen our three benchmark scenarios indicated by 
the sohd green hues in the right panel of Fig. [H which in principle correspond to the three 
basic situations that are possible after the next generation of experiments!^ 



A: negative signal 

In this situation, no signal of Ouf^P has been measured by GERDA. No matter if we 
have additional information on the neutrino mass by kinematic measurements or by 
cosmology, we do not know if neutrinos are Dirac or Majorana particles. Correspond- 
ingly, we cannot use |mee| to constrain and/or S, or vice versa. Concerning new 
physics, we can only give an upper bound on the size of the corresponding amplitude. 

B: consistent positive signal 

Here, GERDA has measured a certain rate of Oz//3/3 and the corresponding range 
obtained for \mf,e\ is consistent with complementary positive signals from KATRIN 
and Planck, as in the "Q'D scenario" of Ref. [SS]. In this case, one can fit all signals 
with the minimal scenario and we expect strong limits on contributions from further 
new physics. Depending on which term dominates, this situation will lead to a strong 
upper bound on |ej| or on the interference term |eicos0j|, cf. discussion in Sec. 15.21 

C: inconsistent positive signal 

Finally, we consider the case where GERDA has measured a non-zero rate of OufSP, 
which does, however, disagree with the values of \mee\ consistent with the measure- 
ments by KATRIN and Planck. In particular, these measurements are not consistent 
with a Oz//3/3 rate as high as measured. Then, we are forced to conclude that Oz//3/3 is 
actually dominated by the contributions from new physics. In such a situation, one 
can estimate or even measure the magnitude of \ei\. 

These three benchmark scenarios will be used in the following to derive the respective 
information on the contributions of point-like new physics to Oi^/3/3. 

Finally, we want to remark that we could, of course, have chosen a different selection 
of experiments, in particular of future and upcoming searches for Oz//3/3 (e.g. EXO [56| [57]. 
Majorana [58], or MOON [59l [60]) or for the kinematical determination of the neutrino 
mass (e.g. MARE [61]). When having a look at the next five years, however, we consider 
our selection of experiments and scenarios to be both, representative and realistic. 

^However, note that there might of course be some subtleties involved, in particular in cases where the 
experimental results point to regions in parameter space which may correspond to more than one of our 
scenarios. We do not take into account such peculiarities, as long as there is no real measured signal which 
indeed falls in exactly such a region. 
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5 Half-lives and numerical analysis 



The full half-life obtained from the Lagrangian in Eq. ([T]) in combination with the standard 
light neutrino exchange can be derived from Refs. [26] and [IT], and is given by 



1=0 



j=4 



Re 




(13) 



The sum consists of three parts, namely the squared terms originating from the amplitudes 
labeled (0,1,2,3), the squared terms originating from the amplitudes labeled (4,5), and 
the interference termsj^ Although all operators have two electrons in the final state, the 
dependence of the nuclear/atomic physics parts on the final state energies are different [IT] , 
which is the reason for the appearance of the three different phase space factors Gi_2,3 [26] . 
The subscript "0" denotes the part coming from light neutrino exchange, while the other 
subscripts label the corresponding operators in Eq. ([1]). Accordingly, one has 



M.Q for % = 0, 
eiMi for i > 0. 



(14) 



Here, the amplitude for the light neutrino exchange is given by the standard expression [T71 
[62], 

2 



Mo = M 



with M 



9 a 



(15) 



and all the other NMEs Mi as well as the phase space factors are specified in Ref. [26] . 

Note that there might be CP- or other phases involved, i.e., the coefficients might 
effectively be complex numbers (where any overall phase will cancel out, but a relative phase 
can remain). Hence, it is necessary to think about how to deal with possible interference 
terms. Sometimes in the literature, such phases are ignored by assuming CP-invariance [53] , 
which then results in relative signs only, which could be plus or minus [19], leading to 
minimal or maximal cancellation of the partial amplitudes. We will discuss this point in 
the following. 

First of all, one should ask the question when interference terms can appear at all. The 
key point is that states of a different helicity will not interfere, as they have a quantum 
number that distinguishes them, a fact that is well known for, e.g., fi~ — )■ e~e~^ scat- 
tering [M]. But how is the situation in a Oi^/3/3-transition? As is usual, we consider only 
a O"*" — 7- O"*" nuclear transition, which is by far the dominant contribution [T7]. Hence, no 
matter which chiralities were chosen in Eq. ([2]), the decay will always have the same initial 
state, namely a nucleus with zero spin and positive parity. Also, the nuclear final state 
will have zero spin. Furthermore, since the s-wave approximation is taken for the final 
state electrons [26], they will have a zero orbital angular momentum. In such a situation, 
the conservation of angular momentum forces the final electron spins to point in opposite 



^Note that, depending on the hehcities of the final state electrons, several of the interference terms in 
Eq. (jl3|) can be zero. This point will be discussed in a moment. 
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directions. However, this does not constrain the hehcities of the final state electrons, since 
the nucleus itself can obtain a non-zero 3-momentum, which allows all electron helicities 
to be producedlll Accordingly, even though we start with the same initial state for all 
transitions, we can have different final states, which will not interfere. 

What we can learn from this discussion is that the decisive point is only the helicities of 
the final state electrons, cf. Eq. ([3]). However, the final state quarks do not matter (exactly 
as the initial state quarks do not matter), since they are "hidden" inside the nucleus, and 
their helicities are never measured. Hence, even though the initial or final state helicities 
in Eq. ([2]) might differ on the quark level, it does not matter for the actual decay, since 
the quarks are only intermediate states. 

This leads us to a question: Which of the operators in Eq. ([1]) can actually have the 
same final state electron helicities as the standard light neutrino exchange process? This 
is easy to answer: Using the symmetries involved, one can see that, in the standard case, 
the electron current is given by Jr from Eq. ([3]), due to the double weak vertex on the 
Majorana neutrino line ^E\. This leads to two left-handed electrons in the final state, as 
to be expected from the left-handedness of the weak interactions]^ 

We can conclude: The only operators that can interfere with the standard process are 
the ones with the coefficients ei, €2, and €3, and even they can only interfere if the current 
jn is involved. Accordingly, one can constrain the following operator coefficients without 
having to care about interferences: e^^^, €3^^, Cg^^, e^^^, and e^^^, where x,y,z G {L,R}. 

Note that interference effects between different modes of Oi^/3/3 have also been discussed 
in Refs. [iQl E]. They both agree with our findings in that the important condition for 
the appearance of interference effects is the equality in the chiral structures of the new 
physics contribution and the light neutrino exchange diagram, and hence the equality of 
the corresponding phase space factors. Furthermore, Ref. [H] points out that the intro- 
duction of new physics contributions could easily lead to relative CP-phases, as we have 
also mentioned. 

Finally, we want to note that interferences will, if they appear at all, only play a role 
in scenario B. The reason is simply that for scenarios A and C, it is perfectly enough to 
consider the contribution from the respective effective operator only. 

These considerations allow us to write down expressions for the half-lives for all three 
benchmark scenarios: 

A: negative signal 

Here, only the contribution from the effective operators can be constrained. The 
expression we need is 

^In other words: The correct hehcities arc enforced by the electrons being emitted in certain directions. 
However, these directions cannot be detected in an experiment hke GERDA. 

^Actually, we arc sloppily setting helicity and chirality equal, which is strictly speaking only justi- 
fied in the highly relativistic limit. Nevertheless, the kinetic energy of both electrons for ^^Ge will be 
2.04 MeV [19], making them relativistic, but not highly relativistic. Although this might lead to some 
error, we have to take into account that the largest error will actually come from the problematic deter- 
mination of the NMEs involved. 
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where i = 1, 2, 3, 4, 5, Gi^2,3 = Gi, and (^4,5 = G2. 

B: consistent positive signal 

In this case, we will also need to consider the contribution coming from the standard 
light neutrino exchange and obtain 

[^vflil = G^\Mof + G,\e,\'\M.f (17) 
(+ 2Gi cos 0i ■ lef^^l lA^ol |M| for z = 1, 2, s) , 

where i = 1,2,3,4,5, and 0j = (p^^^ is the relative phase between the two different 
contributions in the relevant cases. Of course, we cannot give an explicit expression 
for (pi as long as the underlying theory is unknown, but e.g. in the case of the exchange 
of new heavy active neutrinos, this would correspond to a new Majorana phase 7 in 
addition to the two phases a and (3 for light neutrinos. 

C: inconsistent positive signal 

Also here, we only need the contribution from the effective operators, which yields 
the same formula as for scenario A, 

[T'y2''U = G.h\'\M.f- (18) 

These scenarios will now be analyzed numerically. 

5.1 Scenario A: Negative signal 

This scenario is the easiest one to analyze. The key point is that, in the case of a negative 
signal, we cannot use any external information on the neutrino parameters. The reason 
is that, even though we could in principle calculate a range for the effective mass \mee\ 
from data obtained by neutrino oscillation experiments in combination with kinetic and/or 
cosmological neutrino mass measurements, lepton number violation would not have been 
experimentally established, and hence, we would not even know if |mee| has any physical 
meaning. Neutrinos could be Dirac particles, still in agreement with all experimental 
results to datej^ 

The only decent statement we can make in such a situation is that the contribution 
to OuPP from effective point-like operators can be at most as large as the upper bound on 
the decay rate measured in the experiment, i.e., 

r. < Tbound- (19) 

Note that this limit is a conservative one in the sense that we allow the contribution from 
the effective operators to be the dominant one, while any possible interference terms or 

^Of course, this will change if we know from any other type of experiment that lepton number violation 
exists and that neutrinos are Majorana fermions. However, such a situation is not very likely to be present 
before data from the GERDA experiment will be available [4^, so we do not consider this case here. 



10 



A 






1 LLz,HRz\ 


1 LHz,HLz\ 


|e4| 




Bound 


1.6 ■ 10-^ 


9.5-10-11 


1.2 ■ 10-9 


7.3 ■ 10-1° 


7.9 ■ 10-1° 


7.1 ■ 10-9 



Table 1: The resulting upper bounds on the magnitudes of the coefficients of the effective 
operators in scenario A, where z = L ot R. 

further contributions are neglected. One could, of course, improve the bounds by assuming 
additional contributions rendering the point-like interactions much less dominant. How- 
ever, the corresponding limit would only hold under this assumption and would hence be 
much more model-dependent. 

Furthermore, Eq. ( fT9l) should be viewed as an estimate rather than as a strict bound at a 
certain confidence level, due to the intrinsic uncertainties associated with the computation 
of the NMEs. We will comment on this point in Sec. 15.41 

In our analysis, we will always assume a single term from Eq. (fl]) with a certain chirality 
structure to be dominanl0, e.g. e^^^ J^JRiL,!!- This will lead to a certain decay rate T^lrl 
that is a function of the corresponding coefficient e^^^ . Applying Eq. f|T9|) will then yield 
an upper bound on e^^^ . 

A similar analysis has been done in Ref. [26] for a 90 % CL. limit obtained by the 
Heidelberg-Moscow experiment {Ty^^ > 1.8 ■ 10^^ years, Refs. [261 [65]). which we have 
also reproduced as a cross check, using the same numbers as in the original referenced 

To adjust the numbers to our scenario A, we have taken the 90 % CL. limit for GERDA 
phase HI in the case of a negative signal, {Ty^^)x > 2.8 ■ 10^'^ years [S3]. The nuclear 

radius for ^^Ge has been estimated as i? ~ 1.2^/76 fm ~ 5.1 fm, and the numerical values 
for the phase space factors have been taken from Ref. [26]. However, it was necessary to 
update the values of the NMEs, since there have been many developments within the last 
years, in particular in what regards the inclusions of short-range correlations and higher 
order nuclear currents [BSj- Taking the calculations from Refs. [S21ES], we can update the 
values from Ref. [26] to b^ 

Mf,n = 6.8 ■ 10-2 and Mgt,n = 1-7 ■ lO-i. (20) 

Note that we also use qa = 1.25, as suggested in Ref. [62]. In addition, all other numerical 
values are taken from Ref. |26j . 

The resulting upper bounds on the magnitudes of the coefficients of the effective op- 
erators are given in Tab. [H where we have indicated the chirality structures only for the 

^"This is what is called "on-axis evaluation" in Ref. j26j . 

"'^"'^Note that we could reproduce all numbers from Ref. [26] except for the bound on jg^^^^^-f^^j^ where 
our value is smaller by a factor of two. 

^^Note that the matrix elements used in Ref. [IS] were based on the calculations from Ref. [ST], up to 
a different sign convention in the definition of the Fermi part. The values from Refs. [621 [66], however, 
were based on the calculations from Ref. |68| , which uses a normalization for the individual Fermi and the 
Gamow- Teller parts that is smaller by a factor of m^/nip, but the full NME is of course of the same order 
of magnitude. 



11 



cases where the bounds depend on them. As to be expected for scenario A, the bounds on 
the coefficients will be stronger than the bounds based on the Heidelberg-Moscow limit by 
about one order of magnitude, since phase III of GERDA would yield a limit on the decay 
rate that is about two orders of magnitude better, and since the squared absolute values 
of the coefficients enter the decay rate. 

Finally, we would like to also derive a bound on the paremeter A'^^^^^ in our RPV-SUSY 
example. Using Eq. ([6]), as well as the bound on ef^^ from Tab. [H it is easy to derive 
the value A'^^^^^ < 0.06 for the sparticle masses assumed, which is a reasonable value as 
compared to other bounds on this parameter |42j . 



5.2 Scenario B: Consistent positive signal 

In case of a positive signal from GERDA, one first needs to check if all the available data 
from all experiments is consistent with the assumption that light neutrino exchange is 
the only mechanism behind Oz//3/3. We will also investigate situations in which only one 
experiment (KATRIN or Planck) can yield complementary information on the neutrino 
mass, which would be important if one of the two ceases to produce reliable results. 

The main technical problem is how to separate the point-like contributions resulting 
from the interactions in Eq. ([1]) from the dominant contribution of light neutrino exchange, 
which is a long-range force [37]. Note that, in principle, there might also be further long- 
range contributions, such as new heavy W^-bosons in addition to the ones in the SM, in 
which case the neutrino line would still be of long range. However, we do not consider such 
contributions here, but instead focus on the point-like ones that can be fully expressed by 
effective field theory, in order to be as model-independent as possible. 

To perform the analysis, we first need to assume a certain measured rate for scenario B. 
Taking the smallest neutrino mass to be 0.3 eV (cf. Sec. H]), we roughly obtain (by vary- 
ing the neutrino oscillation parameters and the Majorana phases) an effective mass |mee| 
between 0.1 eV and 0.3 eV, which translates into a half-life between 2.8 • 10^^ years and 
2.6 ■ 10^^ years for an NME of 4.0. For definiteness, we assume a certain value within this 
range to be observed, let us say 

T°/f^) 5.0 -lO^^ years. (21) 
This is equivalent to a decay rate of 

ln2 



obs 



1.4 ■ 10"^Vyear. (22) 



A/2 ). 



B 

This measurement will, of course, have a certain error, the estimate of which can be found 
in Ref. [29] to be 

^ 23 %, (23) 

^ obs 

which leads to an absolute error on the decay rate of Uobs = 0.32 ■ 10~^^/year. Finally, 
comparing Eqs. ([7]), ( !T3|) . and (fT5|) . one observes that G = Gi\n2/ml. This enables the 
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Parameter 


Best-fit value & lo" range 


Reference 




0.32 ±0.016 


[28] 


■^13 


0.00 ±0.0076 


m 




(7.6 ±0.19) ■ 10^5 eV^ 


m 


1 ATTig-]^ liior. 


(2.5 ±0.09) ■ 10-^ eV^ 


[28] 


I^^31 |inv. 


(2.3 ±0.10) ■ 10^^ eV^ 


[28] 



Table 2: The best-fit values and (symmetrized) la ranges for the neutrino oscillation 
parameters used in the analysis. Note that, in the inverted hierarchy case, |mee| is actually 
proportional to iy|Am|]^|inv. and depends strongly on 6*12, whereas there is only a very weak 
dependence on 6*13 05]. Accordingly, the former two observables would require a precise 
experimental determination in order to maximize the use of our analysis. 

calculation of the parameters to be used in the x^-function in Eq. flTT]) . namely, from 
Eq. m, 

Imeelobs ^ 0.23 eV, a(|mee|) ^ 0.026 eV, (24) 

where a representative nominal value of [A^^^'l = 4.0 has been used for the NME [291 162] . 
Consistent values for the observables measured by KATRIN and Planck can be found 
in Ref. [SS] for a smallest neutrino mass of 0.3 eV: (m^)obs = (0.30 eV)^ and Sobs = 
0.91 eV, respectively. Realistic values for the corresponding errors are given by cr(mp = 
0.025 eV^ [301 El] and cr(S) = 0.05 eV [IE]. For definiteness, we express all neutrino masses 
by the lightest neutrino mass m using the formulas that hold for an inverted mass ordering. 
However, as we are with m = 0.3 eV in the quasi-degenerate region of the neutrino masses, 
it will anyway not make much of a difference which ordering is used [15]. Finally, the 
values and uncertainties of the neutrino oscillation parameters we use as input are given 
in Tab. [2 

Now we are ready to investigate the x^-function, Eq. (ITT]) , which will only be a function 
of the absolute neutrino mass scale (or, equivalently, the smallest neutrino mass m), the 
two Majorana phases a and /3 (where the dependence on the latter phase is rather weak), 
and the accuracy F to which we know the NME. As already mentioned, we investigate nine 
cases, namely the one where both KATRIN and Planck yield a result as well as the two 
scenarios in which only one of these experiments yields a useful measurement, and each of 
these three cases will be investigated for perfect knowledge of the NME (F = 1) as well as 
for different degrees of imperfect knowledge {F = 2, 3). 

Following the procedure from Ref. [H], we minimize the x^-function with respect to a, 
P, and ^. The results are the functions of m shown in Fig. [21 where we have plotted the 
behavior of the x^-function for the case where only GERDA yields a result (pessimistic), 
as well as for the cases where also KATRIN and/or Planck yield results. Note that a 
measurement of 0^(3(3 only (upper-left panel) does not yield an unambiguous result for the 
neutrino mass scale. This is natural because of the intrinsic uncertainties involved, i.e., 
bad knowledge of the NME involved can, in certain regions of the parameter space, be 
compensated by a variation of the Majorana phase a (and /3, although this dependence 
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Figure 2: The ranges for the smallest neutrino mass ms for different cases. 



is subdominant). Hence, the flat prior for the exact value of the NME corresponds to a 
perfectly flat region in the x^-function0 This is different in the case when all experiments 
yield information (upper- right panel): Then, KATRIN as well as Planck will point to a 
certain neutrino mass scale, and hence, there will be no flat region left in the x^-plot? which 
means that we actually have a measurement of the neutrino mass, and we can determine 
the corresponding best-fit point as well as the la region. The qualitative behavior does 
not change when we have only one experiment in addition to the Oi^/3/3- signal: If there is 
only the kinematical determination (lower-left panel), the x^-function becomes broader, 
but it is nevertheless possible to derive a consistent range for the neutrino mass. The case 

^■^Note that wc display the full x^-function, rather than only Ax^, since our assumption of consistent 
results always allows for (at least) one perfect fit point. 
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F 


KATRIN & Planck 


KATRIN only 


Planck only 


1,2,3 


O.SOlHl^ eV 


OQ4-U.U4U y 


O.SOlHi^ eV 



Table 3: The reconstructed values and la ranges for the smallest neutrino mass. 

of a cosmological observation being the only information in addition to Ou/SP (lower-right 
panel) looks practically identical to the case where all measurements fit together. This 
is simply a reflection of the small uncertainty assigned to the observation of the sum E 
of neutrino masses. However, note that there could be systematic errors involved in the 
cosmological data that we are not aware of, which could even lead to wrong conclusions 
about the actual value of the neutrino mass ^5], which is why this scenario might be a 
little too optimistic in case there is no additional information by KATRIN. 

We have taken into account different values of the degree of knowledge of the NME 
[F = 1,2,3) in all plots. However, for the la region, this only plays a role for the (too) 
pessimistic scenario of having a positive signal from Oz//3/3 only, even though the neutrino 
mass would be in a range where KATRIN and Planck could both determine it. We will 
not consider this case further. 

Using the x^-functions obtained, one can derive the best-fit values and la errors for 
m in all cases, shown in Tab. |3J As to be expected, we always obtain our assumed value 
of m = 0.3 eV as best-fit point, which is a reflection of the fact that we consider a sit- 
uation in which all experiments yield consistent results. Note that, for the la-regions, 
the uncertainty in the NME does not play a major role, since we have enough comple- 
mentary information from the kinematical and from the cosmological determination of the 
neutrino mass. Furthermore, note that, although it seems that the KATRIN experiment 
does not yield much more information compared to using the cosmological data alone, it 
nevertheless cannot be overemphasized that experiments based on kinematics are the only 
model-independent probe for the absolute neutrino mass scale that we have. Depending on 
the outcome of the cosmological observations, it might be necessary to dismiss the data on 
S completely, in case that inconsistencies in our cosmological model would arise. 

Now we are finally in the position to derive bounds on the different |ej|'s. Let us 
shortly summarize what we have done: First of all, we started with Eq. (1T7|1 . where we 
have neglected the interference terms, as discussed at the beginning of Sec. [5], 

[^vf 1b.' = Gi lA^ol' + G,|e,p \Mif . (25) 

Next, we have assumed a certain observed rate, Eq. (1221) . with an error given by a^hs in 
Eq. 0231) . Our goal is to derive bounds on the |ej|'s, so we should re- arrange Eq. (12^ to 
obtain 

GM\' \M.f = [T^f IbJ - Gi lA^ol^ = ^^^^i^- (26) 

Since, as we have shown, the measurement of the neutrino mass scale is consistent with 
that of GERDA, the whole decay rate could be attributed to the standard mechanism of 
light neutrino exchange, causing the right-hand side of Eq. f l26|) to vanish. However, the 
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B 


(!) 


k2| (!) 


1 LLz,HHz\ /t\ 

h \ (!) 


1 LHz,HLz\ /t\ 

h \ (!) 


|e4| 




Bound 


5.8 ■ 10^^ 


3.5- 10^1° 


4.5 ■ 10-9 


2.7- 10-9 


2.9 ■ 10-9 


2.6 ■ 10-s 



Table 4: The resulting upper bounds on the magnitudes of the coefficients of the effective 
operators in scenario B, where z = L or R. Here, "(!)" denotes that there is a possibility 
of interference if the current is involved in the operator. 



effective operators could still contribute within the error of the measurement, giving the 
estimate 



Gi\e,\^\M^^ < 



O'obs 

M2"' 



Finally, Eq. ( 127|) can be rewritten to yield a bound on |ej|: 

, , 1 



0"obs 

a In 2' 



(27) 



(2J 



The numerical results can be found in Tab. |H One observes that the constraints for this 
scenario appear to be slightly worse than the ones from scenario A (cf. Sec. 15.11) . However, 
we can see in Fig. [1] that, for scenario A, we would actually need a successful phase HI 
of GERDA, while a situation like scenario B could even be reached within phase I, the 
sensitivity of which is worse by a factor of roughly 100. Since the decisive quantity for 
the bounds is the error rather than the decay rate itself, cf. Eq. ( l28l) . the bounds obtained 
would be strong even for phase I. However, this is only possible if there is no inconsistency 
in the experimental results, cf. Sec. 15.31 

Using the example from Eq. ([6]) again, we now arrive at A'^^^^ < 0.12, which is slightly 
worse than the value from Sec. 15.11 due to the difference on the respective bounds (cf. 
Tabs. □ and ij). 

One last important remark is that some of our results, i.e., those marked with "(!)" 
in Tab. HI could be altered by the presence of an interference term. However, as already 
discussed, these terms will only play a role when two conditions are simultaneously fulfilled: 

1. the effective operator contains the electron current (otherwise there can be no 
interference in the highly relativistic limit), and 

2. the relative phase (j) between the standard operator and the new physics contribution 
is not close to 7r/2. 

The ffist condition is fulfilled for only a handful operators [cf. Eqs. ([T]) and ([3])], and even 
then there is still the possibility of cancellation due to the value of the relative phase 
(second condition). Nevertheless, we also give the bounds for such a situation. In those 
cases where interferences might be relevant, we expect that, unless cos (pi is very small in 
Eq. ( flTl) . the last term will dominate over the second one. Then, once again assuming 
that the standard mechanism of light neutrino exchange is the main contribution, one can 
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B (interf.) 


COS 01 


1 xyH 1 1 
e2 COS 02 


1 LiLiH .tittti 1 1 
63 cos 03 


1 LHH,RLH , 1 
63 COS 03 


Bound 


1.4 ■ 10^8 


8.4 ■ 10^11 


1.1 ■ 10^9 


6.5 ■ 10^1° 



Table 5: The resulting upper bounds on the combination of parameters |ef^ cos0i| in 
scenario B, assuming that the interference terms dominate. Here, x,y = L or R. 

calculate an approximate value for \Aio\ and then derive bounds on the combination of 
parameters |e^^^cos0j| as 

lef^"^ cos 0,1 < ^0.245,, 2 = 1,2,3, (29) 

V obs 

where Bi is the left-hand side of Eq. f l28|l . the values of which are given in Tab. HI The 
resulting bounds on |e^^'^cos0j| are shown in Tab. [51 These bounds are seemingly stronger 
by a factor of roughly four than the ones without interferences, but one has to keep in 
mind that the cosines involved may very well be of the order of 0.1, which would, in turn, 
weaken the bounds on the |e^^'^|'s. 

Here, the RPV-SUSY parameter can only be constrained in connection with a CP-phase 
0A [ll]- The corresponding bound from Tab. [5] translates into |A']^^]^ cos 0a| < 0.028, which 
could turn out to be particularly strong in case the value of the CP-phase is predicted in 
a certain model. 

5.3 Scenario C: Inconsistent positive signal 

The final case to consider is a positive signal in GERDA whose rate cannot be described 
by neutrino physics only. As an example, we take a measured decay rate of {T^l2^)c = 
1.26 ■ 10^^ years (well below the GERDA phase H limit), which would translate into an 
effective mass of \mee\ = 0.045 eV for an NME of 4.0. In case of a negative signal in 
KATRIN and Planck [and under the implicit assumption that there are no unknowns in 
cosmology which would affect S in a way that it is not given by the expression in Eq. ( ITOj) 
anymore], we would practically exclude inverted mass ordering, since the minimum value of 
S in this case is Sj^j^^ = 0.1 eV, which is just excluded by a negative signal from cosmology. 

In this case, the contribution of the light neutrino exchange to the decay rate would be 
smaller than the contribution from the effective operators by some factor. Then, under the 
assumption that the neutrino contribution can be completely neglected, one can actually 
derive a value (and, in principle, allowed ranges), which serves as estimate for the coeffi- 
cients €{. Note that this would be in particular valid in case that the parameters are such 
that |mee| is practically zero [15]. In Fig. [3l the minimized x^-function as a function of the 
smallest neutrino mass is shown for normal and inverted mass ordering. We can see that, 
although we could in principle explain such a high Oi^/3/3 rate for a certain range of the 
neutrino mass, this range does not fit with the additional information obtained by Planck 
which strongly disfavors the mass ranges required by GERDA. Accordingly, the minima of 
the full x^-functions will be at least around 4, indicating that the results are inconsistent. 
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Figure 3: The minimized x^-function as a function of the smallest neutrino mass, which is 
mi for normal mass ordering and for inverted mass ordering. 



We will now use Eq. f lTSj) to estimate the magnitudes of the e^'s. The Oz//3/3 decay rate 
Fobs, cf. Eq. fl22|) . will translate into estimates of the sizes of the e^'s, 




Qlobs ^ -n-TT J ^ . ^ - (30) 



Using the same formulas as in Sec. I5.H we obtain estimates in the case of scenario C shown 
in Tab. [61 

We conclude that one could, in principle, constrain the coefficients of the effective 
operators very well by simply assuming that the contributions from neutrino physics are 
completely subdominant. As to be expected, values slightly smaller than the bounds 
obtained in Ref. |26] could be measured. 



18 



c 




1^2! 


1 LLz.HHzi 


1 LHz,HLzi 


|e4| 




Estimate 


2.3 ■ 10-^ 


1.4-10-1° 


1.8 ■ 10-9 


1.1 ■ 10-9 


1.2 ■ 10-9 


1.1 ■ 10"^ 



Table 6: The resulting estimates for the magnitudes of the coefficients of the effective 
operators in scenario C, where z = L ot R. 

Accordingly, one can also derive the magnitude of the RPV-SUSY parameter in this 
scenario, which turns out from Tab. [6] to be A'^^^^ ^ 0.074. In principle, this method could 
be used to measure the amount of i?-parity violation present. 

5.4 On the importance of the knowledge of the nuclear physics 

There is an important caveat in our argumentation: In the above analyses, we have (in 
connection with the effective operators) always used the values for the NMEs derived from 
Refs. [62l [66] (which update the ones given in Ref. [26]). However, in reality these values 
will be uncertain. Different working groups might obtain different numbers, which could 
then be compared to the experimental results. A more advanced analysis that takes this 
into account has been performed in our scenarios B and C for the light neutrino exchange 
diagram. However, for the effective operators we do not want to guess some values of the 
NMEs that may be more or less likely, as long as no further calculations exist. One still 
has to keep in mind that exactly these computations might become very important in case 
that, after the next generation of experiments will have been completed, we would indeed 
end up in a situation that is similar to one of our scenarios. 

The goal of this paper was to demonstrate how strong the bounds on point-like new 
physics can become in the case of a positive or negative future signal of Oz//3/3. This has been 
achieved by the estimates obtained in Sees. 15.11 15. 2[ and 15.31 However, this does not cure 
the problem of the relatively poor knowledge of the nuclear physics involved. Although 
this might not matter too much in a phenomenological analysis like the one presented 
here, it will be crucial for the true evaluation of future data - a point that cannot be 
overemphasized. Hence, our results can very well serve as a guideline for how existing 
bounds on the effective operators will improve for different realistic situations. However, 
the numbers obtained suffer from the uncertainties involved in the computations of the 
nuclear matrix elements and can, accordingly, not be treated as strict statistical bounds. 

6 Conclusions 

We have investigated point-like operators contributing to neutrino-less double beta decay 
using the formalism of effective field theory. In order to focus on situations that can be 
tested in the near future, we have defined three scenarios. A, B, and C, which correspond 
to a negative double beta signal, a positive double beta signal that is consistent with com- 
plementary experimental results, and a positive double beta signal that is inconsistent with 
complementary experimental results, respectively. For each scenario, we have determined 
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the bounds on (or estimates for) the strengths of each point-hke operator that could con- 
tribute to the decay. Our scenario A is merely an update of bounds that have already been 
present in the literature, while our scenarios B and C have, to our knowledge, not been 
investigated before in this context. In case neutrino-less double beta decay is not detected 
at all, the current bounds on the e^'s can be improved by roughly one order of magnitude. 
The most interesting case arises if a positive signal is consistent with all complementary 
results: Then, the possible room for further new physics contributions to double beta decay 
on top of the standard light neutrino exchange is essentially determined by the error in 
the measured rate, which is much smaller than the rate itself. In case of an inconsistent 
positive signal, one can even give estimates of the coefficients €{. Our results do, however, 
suffer from the lack of knowledge of the underlying nuclear physics, which is a well-known 
problem in studies of neutrino-less double beta decay. Accordingly, the bounds that we 
have derived will be made much more robust once the nuclear matrix elements are known 
to a better precision. In that case, our bounds would depend directly on the experimental 
accuracy and could be significantly improved by future attempts to measure the rate of 
neutrino-less double beta decay more precisely. 
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